Feasibility of proton dosimetry overriding planning CT with daily CBCT elaborated through generative artificial intelligence tools.

Radiotherapy commonly utilizes cone beam computed tomography (CBCT) for patient positioning and treatment monitoring. CBCT is deemed to be secure for patients, making it suitable for the delivery of fractional doses. However, limitations such as a narrow field of view, beam hardening, scattered radiation artifacts, and variability in pixel intensity hinder the direct use of raw CBCT for dose recalculation during treatment. To address this issue, reliable correction techniques are necessary to remove artifacts and remap pixel intensity into Hounsfield Units (HU) values. This study proposes a deep-learning framework for calibrating CBCT images acquired with narrow field of view (FOV) systems and demonstrates its potential use in proton treatment planning updates. Cycle-consistent generative adversarial networks (cGAN) processes raw CBCT to reduce scatter and remap HU. Monte Carlo simulation is used to generate CBCT scans, enabling the possibility to focus solely on the algorithm's ability to reduce artifacts and cupping effects without considering intra-patient longitudinal variability and producing a fair comparison between planning CT (pCT) and calibrated CBCT dosimetry. To showcase the viability of the approach using real-world data, experiments were also conducted using real CBCT. Tests were performed on a publicly available dataset of 40 patients who received ablative radiation therapy for pancreatic cancer. The simulated CBCT calibration led to a difference in proton dosimetry of less than 2%, compared to the planning CT. The potential toxicity effect on the organs at risk decreased from about 50% (uncalibrated) up the 2% (calibrated). The gamma pass rate at 3%/2 mm produced an improvement of about 37% in replicating the prescribed dose before and after calibration (53.78% vs 90.26%). Real data also confirmed this with slightly inferior performances for the same criteria (65.36% vs 87.20%). These results may confirm that generative artificial intelligence brings the use of narrow FOV CBCT scans incrementally closer to clinical translation in proton therapy planning updates.

Retrospective reconstruction of four-dimensional magnetic resonance from interleaved cine imaging - A comparative study with four-dimensional computed tomography in the lung.

Background and purpose: Imaging of respiration-induced anatomical changes is essential to ensure high accuracy in radiotherapy of lung cancer. We expanded here on methods for retrospective reconstruction of time-resolved volumetric magnetic resonance (4DMR) of the thoracic region and benchmarked the results against 4D computed tomography (4DCT). Materials and method: MR data of six lung cancer patients were collected by interleaving cine-navigator images with 2D data frame images, acquired across the thorax. The data frame images have been stacked in volumes based on a similarity metric that considers the anatomical deformation of lungs, while addressing ambiguities in respiratory phase detection and interpolation of missing data. The resulting images were validated against cine-navigator images and compared to paired 4DCTs in terms of amplitude and period of motion, assessing differences in internal target volume (ITV) margin definition. Results: 4DMR-based motion amplitude was on average within 1.8 mm of that measured in the corresponding 2D cine-navigator images. In our dataset, the 4DCT motion and the 4DMR median amplitude were always within 3.8 mm. The median period was generally close to CT references, although deviations up to 24 % have been observed. These changes were reflected in the ITV, which was generally larger for MRI than for 4DCT (up to 39.7 %). Conclusions: The proposed algorithm for retrospective reconstruction of time-resolved volumetric MR provided quality anatomical images with high temporal resolution for motion modelling and treatment planning. The potential for imaging organ motion variability makes 4DMR a valuable complement to standard 4DCT imaging.

Extension of the cone-beam CT field-of-view using two complementary short scans.

BACKGROUND: Robotic C-arm cone-beam computed tomography (CBCT) scanners provide fast in-room imaging in radiotherapy. Their mobility extends beyond performing a gantry rotation, but they might encounter obstructions to their motion which limit the gantry angle range. The axial field-of-view (FOV) of a reconstructed CBCT image depends on the acquisition geometry. When imaging a large anatomical location, such as the thorax, abdomen, or pelvis, a centered cone beam might be insufficient to acquire untruncated projection images. Some CBCT scanners can laterally displace their detector and collimate the beam to increase the FOV, but the gantry must then perform a 360° rotation to provide complete data for reconstruction. PURPOSE: To extend the FOV of a CBCT image with a single short scan (gantry angle range of 180 ∘ + $180^{\circ}+$ fan angle) using two complementary short scans. METHODS: We defined an acquisition protocol using two short scans during which the source follows the same trajectory and where the detector has equal and opposite tilt and/or offset between the two scans, which we refer to as complementary scans. We created virtual acquisitions using a Monte Carlo simulator on a digital anthropomorphic phantom and on a computed tomography (CT) scan of a patient abdomen. For our proposed method, each simulation produced two complementary sets of projections, which were weighted for redundancies and used to reconstruct one CBCT image. We compared the resulting images to the ground truth phantoms and simulations of conventional scans. RESULTS: Reconstruction artifacts were slightly more prominent in the complementary scans w.r.t. a complete scan with untruncated projections but matched those in a single short scan without truncation. When analyzing reconstructed scans from simulated projections with scatter and corrected with prior CT information, we found a global agreement between complementary and conventional scan approaches. CONCLUSIONS: When dealing with a limited range of motion of the gantry of a CBCT scanner, two complementary short scans are a technically valid alternative to a full 360° scan with equal FOV. This approach enables FOV extension without collisions or hardware upgrades.

Multi-camera optical tracking and fringe pattern analysis for eye surface profilometry in ocular proton therapy.

Background and purpose: An optical tracking system for high-precision measurement of eye position and orientation during proton irradiation of intraocular tumors was designed. The system performed three-dimensional (3D) topography of the anterior eye segment using fringe pattern analysis based on Fourier Transform Method (FTM). Materials and methods: The system consisted of four optical cameras and two projectors. The design and modifications to the FTM pipeline were optimized for the realization of a reliable measurement system. Of note, phase-to-physical coordinate mapping was achieved through the combination of stereo triangulation and fringe pattern analysis. A comprehensive pre-clinical validation was carried out. Then, the system was set to acquire the eye surface of patients undergoing proton therapy. Topographies of the eye were compared to manual contouring on MRI. Results: Pre-clinical results demonstrated that 3D topography could achieve sub-millimetric accuracy (median:0.58 mm) and precision (RMSE:0.61 mm) in the clinical setup. The absolute median discrepancy between MRI and FTM-based anterior eye segment surface reconstruction was 0.43 mm (IQR:0.65 mm). Conclusions: The system complied with the requirement of precision and accuracy for image guidance in ocular proton therapy radiation and is expected to be clinically tested soon to evaluate its performance against the current standard.

The Role of Carbon Ion Therapy in the Changing Oncology Landscape-A Narrative Review of the Literature and the Decade of Carbon Ion Experience at the Italian National Center for Oncological Hadrontherapy.

BACKGROUND: Currently, 13 Asian and European facilities deliver carbon ion radiotherapy (CIRT) for preclinical and clinical activity, and, to date, 55 clinical studies including CIRT for adult and paediatric solid neoplasms have been registered. The National Center for Oncological Hadrontherapy (CNAO) is the only Italian facility able to accelerate both protons and carbon ions for oncological treatment and research. METHODS: To summarise and critically evaluate state-of-the-art knowledge on the application of carbon ion radiotherapy in oncological settings, the authors conducted a literature search till December 2022 in the following electronic databases: PubMed, Web of Science, MEDLINE, Google Scholar, and Cochrane. The results of 68 studies are reported using a narrative approach, highlighting CNAO's clinical activity over the last 10 years of CIRT. RESULTS: The ballistic and radiobiological hallmarks of CIRT make it an effective option in several rare, radioresistant, and difficult-to-treat tumours. CNAO has made a significant contribution to the advancement of knowledge on CIRT delivery in selected tumour types. CONCLUSIONS: After an initial ramp-up period, CNAO has progressively honed its clinical, technical, and dosimetric skills. Growing engagement with national and international networks and research groups for complex cancers has led to increasingly targeted patient selection for CIRT and lowered barriers to facility access.

A new system for in-room image guidance in particle therapy at CNAO.

This paper describes the design, installation, and commissioning of an in-room imaging device developed at the Centro Nazionale di Adroterapia Oncologica (CNAO, Pavia, Italy). The system is an upgraded version of the one previously installed in 2014, and its design accounted for the experience gained in a decade of clinical practice of patient setup verification and correction through robotic-supported, off-isocenter in-room image guidance. The system's basic feature consists of image-based setup correction through 2D/3D and 3D/3D registration through a dedicated HW/SW platform. The major update with respect to the device already under clinical usage resides in the implementation of a functionality for extending the field of view of the reconstructed Cone Beam CT (CBCT) volume, along with improved overall safety and functional optimization. We report here details on the procedures implemented for system calibration under all imaging modalities and the results of the technical and preclinical commissioning of the device performed on two different phantoms. In the technical commissioning, specific attention was given to the assessment of the accuracy with which the six-degrees-of-freedom correction vector computed at the off-isocenter imaging position was propagated to the planned isocentric irradiation geometry. During the preclinical commissioning, the entire clinical-like procedure for detecting and correcting imposed, known setup deviation was tested on an anthropomorphic radioequivalent phantom. Results showed system performance within the sub-millimeter and sub-degree range according to project specifications under each imaging modality, making it ready for clinical application.

Optimizing b-values schemes for diffusion MRI of the brain with segmented Intravoxel Incoherent Motion (IVIM) model.

PURPOSE: To define an optimal set of b-values for accurate derivation of diffusion MRI parameters in the brain with segmented Intravoxel Incoherent Motion (IVIM) model. METHODS: Simulations of diffusion signals were performed to define an optimal set of b-values targeting different perfusion regimes, by relying on an optimization procedure which minimizes the total relative error on estimated IVIM parameters computed with a segmented fitting procedure. Then, the optimal b-values set was acquired in vivo on healthy subjects and skull base chordoma patients to compare the optimized protocol with a clinical one. RESULTS: The total relative error on simulations decreased of about 40% when adopting the optimal set of 13 b-values (0 10 20 40 50 60 200 300 400 1200 1300 1400 1500 s/mm2 ), showing significant differences and increased precision on D and f estimates with respect to simulations with a non-optimized b-values set. Similarly, in vivo acquisitions demonstrated a dependency of IVIM parameters on the b-values array, with differences between the optimal set of b-values and a clinical non-optimized acquisition. IVIM parameters were compatible to literature values, with D (0.679/0.701 [0.022/0.008] ·10-3 mm2 /s), f (5.49/5.80 [0.70/1.14] %), and D* (8.25/7.67 [0.92/0.83] ·10-3 mm2 /s) median [interquartile range] estimates for white matter/gray matter in volunteers and D (0.709/0.715/1.06 [0.035/0.023/0.271] ·10-3 mm2 /s), f (7.08/7.84/21.54 [1.20/1.06/6.05] %), and D* (10.85/11.84/2.32 [1.38/2.32/4.94] ·10-3 mm2 /s) for white matter/gray matter/Gross Tumor Volume in patients with skull-base chordoma tumor. CONCLUSIONS: The definition of an optimal b-values set can improve the estimation of quantitative IVIM parameters. This allows setting up an optimized approach that can be adopted for IVIM studies in the brain.

Synthetic CT in Carbon Ion Radiotherapy of the Abdominal Site.

The generation of synthetic CT for carbon ion radiotherapy (CIRT) applications is challenging, since high accuracy is required in treatment planning and delivery, especially in an anatomical site as complex as the abdomen. Thirty-nine abdominal MRI-CT volume pairs were collected and a three-channel cGAN (accounting for air, bones, soft tissues) was used to generate sCTs. The network was tested on five held-out MRI volumes for two scenarios: (i) a CT-based segmentation of the MRI channels, to assess the quality of sCTs and (ii) an MRI manual segmentation, to simulate an MRI-only treatment scenario. The sCTs were evaluated by means of similarity metrics (e.g., mean absolute error, MAE) and geometrical criteria (e.g., dice coefficient). Recalculated CIRT plans were evaluated through dose volume histogram, gamma analysis and range shift analysis. The CT-based test set presented optimal MAE on bones (86.03 ± 10.76 HU), soft tissues (55.39 ± 3.41 HU) and air (54.42 ± 11.48 HU). Higher values were obtained from the MRI-only test set (MAEBONE = 154.87 ± 22.90 HU). The global gamma pass rate reached 94.88 ± 4.9% with 3%/3 mm, while the range shift reached a median (IQR) of 0.98 (3.64) mm. The three-channel cGAN can generate acceptable abdominal sCTs and allow for CIRT dose recalculations comparable to the clinical plans.

Microstructural parameters from DW-MRI for tumour characterization and local recurrence prediction in particle therapy of skull-base chordoma.

BACKGROUND: Quantitative imaging such as Diffusion-Weighted MRI (DW-MRI) can be exploited to non-invasively derive patient-specific tumor microstructure information for tumor characterization and local recurrence risk prediction in radiotherapy. PURPOSE: To characterize tumor microstructure according to proliferative capacity and predict local recurrence through microstructural markers derived from pre-treatment conventional DW-MRI, in skull-base chordoma (SBC) patients treated with proton (PT) and carbon ion (CIRT) radiotherapy. METHODS: Forty-eight patients affected by SBC, who underwent conventional DW-MRI before treatment and were enrolled for CIRT (n = 25) or PT (n = 23), were retrospectively selected. Clinically verified local recurrence information (LR) and histological information (Ki-67, proliferation index) were collected. Apparent diffusion coefficient (ADC) maps were calculated from pre-treatment DW-MRI and, from these, a set of microstructural parameters (cellular radius R, volume fraction vf, diffusion D) were derived by applying a fine-tuning procedure to a framework employing Monte Carlo simulations on synthetic cell substrates. In addition, apparent cellularity (ρapp ) was estimated from vf and R for an easier clinical interpretation. Histogram-based metrics (mean, median, variance, entropy) from estimated parameters were considered to investigate differences (Mann-Whitney U-test, α = 0.05) in estimated tumor microstructure in SBCs characterized by low or high cell proliferation (Ki-67). Recurrence-free survival analyses were also performed to assess the ability of the microstructural parameters to stratify patients according to the risk of local recurrence (Kaplan-Meier curves, log-rank test α = 0.05). RESULTS: Refined microstructural markers revealed optimal capabilities in discriminating patients according to cell proliferation, achieving best results with mean values (p-values were 0.0383, 0.0284, 0.0284, 0.0468, and 0.0088 for ADC, R, vf, D, and ρapp, respectively). Recurrence-free survival analyses showed significant differences between populations at high and low risk of local recurrence as stratified by entropy values of estimated microstructural parameters (p = 0.0110). CONCLUSION: Patient-specific microstructural information was non-invasively derived providing potentially useful tools for SBC treatment personalization and optimization in particle therapy.

Impact of secondary particles on the magnetic field generated by a proton pencil beam: a finite-element analysis based on Geant4-DNA simulations.

PURPOSE: To investigate the static magnetic field generated by a proton pencil beam as a candidate for range verification by means of Monte Carlo simulations, thereby improving upon existing analytical calculations. We focus on the impact of statistical current fluctuations and secondary protons and electrons. METHODS: We considered a pulsed beam (10 µ ${\umu}$ s pulse duration) during the duty cycle with a peak beam current of 0.2 µ $\umu$ A and an initial energy of 100 MeV. We ran Geant4-DNA Monte Carlo simulations of a proton pencil beam in water and extracted independent particle phase spaces. We calculated longitudinal and radial current density of protons and electrons, serving as an input for a magnetic field estimation based on a finite element analysis in a cylindrical geometry. We made sure to allow for non-solenoidal current densities as is the case of a stopping proton beam. RESULTS: The rising proton charge density toward the range is not perturbed by energy straggling and only lowered through nuclear reactions by up to 15%, leading to an approximately constant longitudinal current. Their relative low density however (at most 0.37 protons/mm3 for the 0.2  µ ${\umu}$ A current and a beam cross-section of 2.5 mm), gives rise to considerable current density fluctuations. The radial proton current resulting from lateral scattering and being two orders of magnitude weaker than the longitudinal current is subject to even stronger fluctuations. Secondary electrons with energies above 10 eV, that far outnumber the primary protons, reduce the primary proton current by only 10% due to their largely isotropic flow. A small fraction of electrons (<1%), undergoing head-on collisions, constitutes the relevant electron current. In the far-field, both contributions to the magnetic field strength (longitudinal and lateral) are independent of the beam spot size. We also find that the nuclear reaction-related losses cause a shift of 1.3 mm to the magnetic field profile relative to the actual range, which is further enlarged to 2.4 mm by the electron current (at a distance of ρ = 50 $\rho =50$  mm away from the central beam axis). For ρ > 45 $\rho >45$  mm, the shift increases linearly. While the current density variations cause significant magnetic field uncertainty close to the central beam axis with a relative standard deviation (RSD) close to 100%, they average out at a distance of 10 cm, where the RSD of the total magnetic field drops below 2%. CONCLUSIONS: With the small influence of the secondary electrons together with the low RSD, our analysis encourages an experimental detection of the magnetic field through sensitive instrumentation, such as optical magnetometry or SQUIDs.

Exploring the role of neutrophil-to-lymphocyte ratio and blood chemistry in head and neck adenoid cystic carcinomas treated with carbon ion radiotherapy.

BACKGROUND AND PURPOSE: In recent years, there is an emerging interest in the prognostic role of chemistry blood biomarkers in oncological patients but their role in adenoid cystic carcinomas (ACCs) is still unknown. This study aims to assess the prognostic significance of baseline neutrophil-to-lymphocyte ratio (NLR) and blood chemistry in a series of head and neck ACC patients treated with carbon ion radiotherapy (CIRT). MATERIAL AND METHODS: We retrospectively retrieved the data of 49 consecutive head and neck ACC patients treated with CIRT. Univariable and multivariable Cox proportional hazard regression (Cox-ph) analyses were performed to look for a potential association of NLR, and other blood biomarker values, with disease-free survival (DFS), Local Control (LC), Metastasis Free Survival (MFS) and overall survival (OS). RESULTS: No significant association between NLR > 2,5 and DFS, LC, MFS and OS was found with univariable analysis although a trend was reported for DFS (Hazard ratio [HR]: 2,10, 95 % CI: 0,85 - 5,08, p-value = 0,11). Patients with hemoglobin (hb) ≤ 14 g/dL showed significantly better DFS, MFS and OS. Multivariable regression Cox-ph analysis for DFS, adjusted for margin status, clinical target volume and Absolute Number of Monocytes, reported the following statistically significant HRs, for both NLR > 2,5 and hb > 14 g/dL respectively: 4,850 (95 % CI = 1,408 - 16,701, p = 0,012) and 3,032 (95 % CI = 1,095 - 8,393, p = 0,033). Moreover, hb > 14 with HR = 3,69 (95 % CI: 1,23 - 11,07, p-value = 0,02), was a negative independent prognostic predictor for MFS. CONCLUSIONS: Pre-treatment NLR and hb values seem to be independent prognostic predictor for clinical outcomes in head and neck ACC patients. If their role will be validated in a larger prospective cohort, they might be worthwhile for a pre-treatment risk stratification in patients treated with CIRT.

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO.

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam's Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

MRI and FUNDUS image fusion for improved ocular biometry in Ocular Proton Therapy.

INTRODUCTION: Ocular biometry in Ocular Proton Therapy (OPT) currently relies on a generic geometrical eye model built by referencing surgically implanted markers. An alternative approach based on image fusion of volumetric Magnetic Resonance Imaging (MRI) and panoramic fundus photography was investigated. MATERIALS AND METHODS: Eighteen non-consecutive uveal melanoma (UM) patients, who consented for an MRI and had their tumour base visible on panoramic fundus photography, were included in this comparative analysis. Through generating digitally-reconstructed projections from MRI images using the Lambert azimuthal equal-area projection, 2D-3D image fusion between fundus photography and an eye model delineated on MRI scans was achieved and allowed for a novel definition of the target base (MRI + FCTV). MRI + FCTV was compared with MRI-only delineation (MRIGTV) and the conventional (EyePlan) target definition (EPCTV). RESULTS: The combined use of fundus photography and MRI to define tumour volumes reduced the average discrepancies by almost 65% with respect to the MRI only tumour definitions when comparing with the conventionally planned EPCTV. With the proposed method, shallow sub-retinal tumour infiltration, otherwise invisible on MRI, can be included in the target volume definition. Moreover, a novel definition of the fovea location improves the accuracy and personalisation of the 3D eye model. CONCLUSION: MRI and fundus image fusion overcomes some of the limitations of ophthalmological MRI for tumour volume definition in OPT. This novel eye tumour modelling method might improve treatment planning personalisation, allowing to better anticipate which patients could benefit from prophylactic treatment protocols for radiation induced maculopathy.

Role of diffusion-weighted MRI in recurrent rectal cancer treated with carbon ion radiotherapy.

Aim: To evaluate the association between pretreatment diffusion-weighted MRI (DW-MRI) and 12-month radiological response in locally recurrent rectal cancer treated with carbon ion radiotherapy. Methods: Histogram analysis was performed on pretreatment DW-MRI for patients re-irradiated with carbon ion radiotherapy for local recurrence of rectal cancer. Results: A total of 17 patients were enrolled in the study. Pretreatment DW-MRI b-value of 1000 s/mm2 (b1000) and apparent diffusion coefficient (ADC) lesion median values for 1-year nonresponders (six patients) and responders (11 patients) demonstrated a median (interquartile of median values) of 62.5 (23.9) and 34.0 (13.0) and 953.0 (277.0) and 942.5 (339.0) µm2/s, respectively. All b1000 histogram features (h-features) and ADC h-kurtosis showed statistically significant differences, whereas only b1000 h-median, b1000 h-interquartile range and ADC h-kurtosis demonstrated remarkable diagnostic accuracy. Conclusion: DW-MRI showed promising results in predicting carbon ion radiotherapy outcome in local recurrence of rectal cancer, particularly with regard to b1000 h-median, b1000 h-interquartile range and ADC h-kurtosis.

Carbon ion radiotherapy is a form of advanced radiotherapy that is especially suitable for radioresistant and/or difficult-to-irradiate tumors. In case of recurrence of rectal cancer after pelvic photon beam radiotherapy, carbon ion radiotherapy may be an option. In this study, the authors looked at the potential role of specific MRI sequences performed before treatment to predict response to carbon ion radiotherapy. If confirmed in a larger prospective cohort, the findings of this study may drive clinical decisions toward a more tumor- and patient-tailored therapeutic approach.

Integrating Biological and Radiological Data in a Structured Repository: a Data Model Applied to the COSMOS Case Study.

Integrating the information coming from biological samples with digital data, such as medical images, has gained prominence with the advent of precision medicine. Research in this field faces an ever-increasing amount of data to manage and, as a consequence, the need to structure these data in a functional and standardized fashion to promote and facilitate cooperation among institutions. Inspired by the Minimum Information About BIobank data Sharing (MIABIS), we propose an extended data model which aims to standardize data collections where both biological and digital samples are involved. In the proposed model, strong emphasis is given to the cause-effect relationships among factors as these are frequently encountered in clinical workflows. To test the data model in a realistic context, we consider the Continuous Observation of SMOking Subjects (COSMOS) dataset as case study, consisting of 10 consecutive years of lung cancer screening and follow-up on more than 5000 subjects. The structure of the COSMOS database, implemented to facilitate the process of data retrieval, is therefore presented along with a description of data that we hope to share in a public repository for lung cancer screening research.

Time-resolved MRI for off-line treatment robustness evaluation in carbon-ion radiotherapy of pancreatic cancer.

PURPOSE: In this study, we investigate the use of magnetic resonance imaging (MRI) for the clinical evaluation of gating treatment robustness in carbon-ion radiotherapy (CIRT) of pancreatic cancer. Indeed, MRI allows radiation-free repeated scans and fast dynamic sequences for time-resolved (TR) imaging (cine-MRI), providing information on inter- and intra-fraction cycle-to-cycle variations of respiratory motion. MRI can therefore support treatment planning and verification, overcoming the limitations of the current clinical standard, that is, four-dimensional computed tomography (4DCT), which describes an "average" breathing cycle neglecting breathing motion variability. METHODS: We integrated a technique to generate a virtual CT (vCT) from 3D MRI with a method for 3D reconstruction from 2D cine-MRI, to produce TR vCTs for dose recalculations. For eight patients, the method allowed evaluating inter-fraction variations at end-exhale and intra-fraction cycle-to-cycle variability within the gating window in terms of tumor displacement and dose to the target and organs at risk. RESULTS: The median inter-fraction tumor motion was in the range 3.33-12.16 mm, but the target coverage was robust (-0.4% median D95% variation). Concerning cycle-to-cycle variations, the gating technique was effective in limiting tumor displacement (1.35 mm median gating motion) and corresponding dose variations (-3.9% median D95% variation). The larger exposure of organs at risk (duodenum and stomach) was caused by inter-fraction motion, whereas intra-fraction cycle-to-cycle dose variations were limited. CONCLUSIONS: This study proposed a method for the generation of TR vCTs from MRI, which enabled an off-line evaluation of gating treatment robustness and suggested its feasibility to support treatment planning of pancreatic tumors in CIRT.

A Dosiomics Analysis Based on Linear Energy Transfer and Biological Dose Maps to Predict Local Recurrence in Sacral Chordomas after Carbon-Ion Radiotherapy.

Carbon Ion Radiotherapy (CIRT) is one of the most promising therapeutic options to reduce Local Recurrence (LR) in Sacral Chordomas (SC). The aim of this work is to compare the performances of survival models fed with dosiomics features and conventional DVH metrics extracted from relative biological effectiveness (RBE)-weighted dose (DRBE) and dose-averaged Linear Energy Transfer (LETd) maps, towards the identification of possible prognostic factors for LR in SC patients treated with CIRT. This retrospective study included 50 patients affected by SC with a focus on patients that presented a relapse in a high-dose region. Survival models were built to predict both LR and High-Dose Local Recurrencies (HD-LR). The models were evaluated through Harrell Concordance Index (C-index) and patients were stratified into high/low-risk groups. Local Recurrence-free Kaplan-Meier curves were estimated and evaluated through log-rank tests. The model with highest performance (median(interquartile-range) C-index of 0.86 (0.22)) was built on features extracted from LETd maps, with DRBE models showing promising but weaker results (C-index of 0.83 (0.21), 0.80 (0.21)). Although the study should be extended to a wider patient population, LETd maps show potential as a prognostic factor for SC HD-LR in CIRT, and dosiomics appears to be the most promising approach against more conventional methods (e.g., DVH-based).

Localization of anatomical changes in patients during proton therapy with in-beam PET monitoring: A voxel-based morphometry approach exploiting Monte Carlo simulations.

PURPOSE: In-beam positron emission tomography (PET) is one of the modalities that can be used for in vivo noninvasive treatment monitoring in proton therapy. Although PET monitoring has been frequently applied for this purpose, there is still no straightforward method to translate the information obtained from the PET images into easy-to-interpret information for clinical personnel. The purpose of this work is to propose a statistical method for analyzing in-beam PET monitoring images that can be used to locate, quantify, and visualize regions with possible morphological changes occurring over the course of treatment. METHODS: We selected a patient treated for squamous cell carcinoma (SCC) with proton therapy, to perform multiple Monte Carlo (MC) simulations of the expected PET signal at the start of treatment, and to study how the PET signal may change along the treatment course due to morphological changes. We performed voxel-wise two-tailed statistical tests of the simulated PET images, resembling the voxel-based morphometry (VBM) method commonly used in neuroimaging data analysis, to locate regions with significant morphological changes and to quantify the change. RESULTS: The VBM resembling method has been successfully applied to the simulated in-beam PET images, despite the fact that such images suffer from image artifacts and limited statistics. Three dimensional probability maps were obtained, that allowed to identify interfractional morphological changes and to visualize them superimposed on the computed tomography (CT) scan. In particular, the characteristic color patterns resulting from the two-tailed statistical tests lend themselves to trigger alarms in case of morphological changes along the course of treatment. CONCLUSIONS: The statistical method presented in this work is a promising method to apply to PET monitoring data to reveal interfractional morphological changes in patients, occurring over the course of treatment. Based on simulated in-beam PET treatment monitoring images, we showed that with our method it was possible to correctly identify the regions that changed. Moreover we could quantify the changes, and visualize them superimposed on the CT scan. The proposed method can possibly help clinical personnel in the replanning procedure in adaptive proton therapy treatments.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors.

PURPOSE: To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario. MATERIALS AND METHODS: The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed. RESULTS: For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D95% and gastro-intestinal OARs D2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs. CONCLUSIONS: The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.